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The induction motor is complex because its dynamic is nonlinear, 
multivariable and highly coupled. In addition, the induction motor (IM) 
controller design can be degraded due to wide ranges operating speed 
specially in low and near-zero speeds and also to variation of interval 
parameters such as the stator resistance change because of stator heating 


during operation. Our contribution is the use of a nonlinear sliding mode 


control using a fuzzy mutual reference adaptive system observer based on 
Keywords: two fuzzy adaptation mechanisms to estimate both rotor speed and stator 
resistance. The proposed work was verified and compared by simulation and 


Fuzzy logic tested in various ranges of speed with motor and regenerating regions of 

Induction motor operation. The results show that the proposed approach has proved good 

Mutual MRAS effectiveness in tracking and observation. 

Sensorless 

Sliding mode This is an open access article under the CC BY-SA license. 
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1. INTRODUCTION 

The industry is often interested in reducing the number of sensors among them the induction motor 
speed sensors. Since the publication of the survey by [1] then [2]-[4], significant interests on sensorless 
sliding mode control of induction motor have been developed in the control research due to the possibility to 
implement this control in industrial applications using new power electronics technology and availability 
cheap and fast computation. Despite the important developments of the sensorless control laws above but 
they suffer from large computation time, complexity, parameters variation, and limited bandwidth control in 
the low and near-zero speed regions estimation. 

In this paper, a nonlinear sliding mode speed and rotor flux control combined with a stable fuzzy 
mutual model reference adaptive system (FM-MRAS) speed and stator resistance estimation for the induction 
motor drive has been developed. By combination the ideas of [5]-[8], the proposed FM-MRAS utilizes two 
models called as reference and adjustable models whose based on the rotor flux errors vectors [9]-[12]. Then, 
these errors are used drive adaptation mechanisms based on the fuzzy logic that generates in parallel the 
estimated speed and stator resistance, this parallel estimation is designed based on the concept of 
hyperstability [6] in order to make the system asymptotically stable. The important advantage of using fuzzy 
logic is robustness against parameter variations. The whole proposed system robustness, stability, 
effectiveness in tracking is tested in simulation. Finally, some comments are presented. 
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2. MATHEMATICAL MODEL OF INDUCTION MOTOR 
The model used is the model whose electrical quantities are all expressed in a fixed frame linked to 


the stator and which will be qualified model(«, 8). The model is given by [13], [14]. 


x=fix)t+g.u 

y = O(x) (1) 
With 

= [Usa Usp]? (2) 

x= lisa isp Yra Prg 2,]7 = [x1 X2 X3 X4 Xs] 


The state vector x is composed of two electrical states (isa isp) and two magnetic states 
(Pra Wrp) and a state representing the rotational speed of the rotor 2,, this speed will be presented by a 
mechanical equation f, gand Lare functions of the variables x .Hence the nonlinear system is of order 5. 
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Where R,is the stator resistance; R,is the rotor resistance; L,is the stator inductance; L,is the rotor 
inductance; L,,is the mutual inductance; T;rotor time constant; ø is the leakage coefficient; pis the number of 
pole pairs; J is moment of inertia; C,is the load torque. 


3. NONLINEAR SLIDING MODE CONTROL 
3.1. Search for a control ensuring a sliding mode 
Definition 4. The control of Figure 1 presents a problem of chattering effects due to switching 


delays. To remedy this problem, we must choose a switching curve called sliding surface S(t). The goal we 
are looking for is to join from the initial state the curve, then force the trajectory to stay on this surface. 


s(t) = (4+ k)" e (4) 


where e is the tracking error that it is equal to the difference between system state space reference and the 
system state space (e = X;e¢ — x), K is the weighting factor that must be strictly positive and n is the degree 
of sliding mode surface [2]. Then, the dynamics of the system will be fixed by the equation of the sliding 
surface, and the origin will be joined with the time constant k. 


S=ke+e=0 


Figure. Chattering effects and switching system on the sliding surface 
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Theorem 4. (Globalstability) assume that there exists a scalar function V (Lyapunov function) of the 
state x, with continuous first order derivatives such that. 


- V(x) = Es 2is positive definite. 

- V(x) = SSis negative definite. 

- V(x) > was lel =>% 

then the equilibrium at the origin S = 0 is globally asymptotically stable [2]. In other words, to follow the 
desired state, it’s enough to make the sliding surface attractive (SS < 0) and invariant (S = 0) [15]. To 


eliminate the uncertain effect of the model and external disturbance, the sliding mode control law is defined 
by adding a command commuting U,to the equivalent command U,,in the form [16]-[18]. 


U = Ueq + Ue (5) 
with 

U; = c sign (S) (6) 
c is a constant and is chosen enough positive to satisfy attractivity and stability conditions. 
3.2. Selection of the control law 

Our goal is to generate a control law [Usa UsgļTusing the sliding control theory. The states 
considered for the induction motor control are: Speed and rotor flux (2,, ¥,-). In the control of the machine, it 


is necessary to ensure the control of the flux then of the speed, when the dynamics of speed is the slowest. 
For convenience of calculation, we will take [19], [20]. 


Y. = Wr = Pha + Wig = x5 + x5 (7) 


Its corresponding derivative gives. 


Ww. = 2X3X3 + 2X4X4 (8) 


According to the (4) the choice of the sliding surfaces of the variables to controlare given in the 
form of a first-order surface system. 


Sı = ky (as — Or ref) +s > Arner) (9) 
So = ka (H, — Brne) +H Bror) 
The successive derivatives of S}andS, will be 
$ = kı (ž5 — Bror) Xe Brne) (10) 
S, = k (Ë, a Pop) + (Ë, = Pep) 
Using the system (3) and (10), the development of calculus gives us. 
. 1 1 : s 
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Such as. 


Xp = XX3 + XOX 
1X) = XX3 + XX4 


S. 2 
X = Xi + x3 
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The system (11) can be rearranged into the following matrix form. 


Sı S F + —Ox 4 OX 3 Usa (12) 
S, F, ÔLnX3 ÖLnX4 || usg 
5 a ae 
with 
Fy = (kı - Ty -42 - ki — pxs(S, + FY) - kô ref — Pr 
2+LmK “i (13) 
F, = Cg Tk)», + LMG $3 + (k2 = 5 aye pat 


To check the Lyaponov’s stability condition, one must have. 

$ = —c sign(S) (14) 
Equalizing (12) and (14) we will have. 

U = —D™tc sign(S) — D71F (15) 


According to (5) the control law components are given by. 


Usg = —D-1F 
sina; (16) 
U, = —D™tc sign(S) 


4. FUZZYMUTUALMRASOBSERVER 
4.1. MRAS based rotor speed estimation 

The principle of speed estimation by MRAS based on rotor flux consists in comparing the outputs of 
two flux estimators obtained by different formulations. The first estimator is based on the stator equations 
defining a model in voltage called the reference model. This model is used to estimate the components flux 
Wray and w,g,from the measurement of voltages and stator currents; being independent of the speed to 
estimate it according to the following [5], [7], [21]. 


apry _ a Pray _ Lr lsa _ a bey 
dt dt A z | ie Usp] Rs [i ahs dt elt a7) 
The second, called adjustable or adaptive model, is described by the rotor equations to estimatethe 
components flux Y,a, and w,g,; Being explicitly linked to speed, he defines a current model is being as. 


Pr Dra ws = Ta m İsa 
Sa siel-b A-6 obiil as) 


Such as @,. is the rotor pulsation with: 0, = p.2,. 

The difference e,,, between the components of the rotor flux delivered at the outputs of the two 
models is used to drive an adaptation mechanism whose objective is to generate the value of the estimated 
speed whitch then reapplied to the input of the adjustable model with. 


Cw, = (Wr py Drar | Pray Drg) (19) 


and 


By = (Kp, +22) ee, (20) 


Figure 2, shows the structure of a classical MRAS speed estimator based on adaptation low generated 
by the adaptation mechanism, which consists of a simple PI-type regulator performing a proportional-integral 
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correction for the induction motor. Although easy to implement this approach, its stability and fast convergence, 
it has some major disadvantages as its dynamic performance sensitivity to parametric variations and the 
accuracy of the estimate is degraded at low and near-zero speeds. The stator resistance variation under operating 
conditions deteriorates the accuracy of the estimate, by modifying at low speed the amplitude and the phase of 
the rotor flux delivered by the reference model. To eliminate those disadvantages, it’s necessary to introduce a 
new method of estimation. This method based on simultaneous and parallel estimation of the rotor speed and 
the stator resistance of the induction motor known as the mutual-MRAS approach. 


[usa al 


[i i | Reference Model 
sarlsg 


' Equation (4.1) 


Ajustable Model Adaptation mechanism 


Equation (4.2) 


Figure 2. Structure of classical MRAS speed estimator based on rotor flux 


In this approach, we introduce a second adaptation mechanism to estimate R,using the same PI-type 
regulator in parallel with the w,adaptaion mechanism, then the error of stator resistance proposed by Popov’s 
hyperstability theoryis represented is being as. 


ers = I sq: (Pray T Pra) F isg- (bg, = Prp,) (21) 
and 


z Kip, 
R, = Ca + zs) er, (22) 
4.2. Proposed observer approach 

In order to improve the performance of the estimation in terms of convergence and robustness, but 
also to avoid the difficulties of setting the PI regulators, we proposed the idea to preserve the basic structure 
of the mutual-MRAS estimator and to use two fuzzy-PI nonlinear controllers forming the adaptation 
mechanisms.This new estimator is of a hybrid nature because, on the one hand, it is based on the synthesis of 
the equations defining the reference model and the adjustable model, and on the other hand it integrates 
control organs derived from artificial intelligence, such us the fuzzy logic. 


4.2.1. Structure of FM-MRAS observer 
The structure of the FM-MRAS observer is illustrated in Figure 3. 


R, (k) 


Fuzzy adaptation 
mechanism for R, 


Equation 
(4.3) 


Equation (4.2) 


Fuzzy adaptation 
mechanism for w, 


Figure 3. Structure of the observer underFM-MRASapproach 
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4.2.2. Synthesis of fuzzy adaptation mechanisms 

The observation of the behavior of the fuzzy mutual-MRAS shows that the significant quantities for 
the simultaneous estimation of the motor speed and the stator resistance are; the error e(k) and its change 
rate Ae(k) for the adaptation mechanism whose are the inputs of the fuzzy correctors.The two fuzzy 
controllers are chosen of identical structure, proportional-integral with incremental output, that is to say that 
their outputs represent the increments of the quantities to estimate w,(k)and R,(k).The structure of the PI- 
Fuzzy adaptation mechanisms is identical for speed and rotor resistance estimation shown in Figure 4.The 
gains Ge, Gaeand Gy, are introduced to normalize the inputs and the output of the estimator controller, they 
are chosen by simulation test to guarantee the belonging of the variables to the associated fuzzy set and thus 
to adjust the sensitivity of the regulator without changing the structure. 


Figure 4. Structure of fuzzy adaptation mechanism 


Normalized inputs are similary determined at time kis being as. 

E(k) = Ge. e(k) (23) 
and 

AE(k) = Gyele(k) — e(k — 1)] (24) 
and the estimated values of the rotor pulsation or the stator resistance are similary deduced by the relation 

u(k) = u(k — 1) + Gau. AU(k) (25) 


or 
u(k) = u(k — 1) + Au(k) (26) 


with u(k ) represents @,(k) or R,(k). 

Fuzzy logic controllers are based on three well known blocs: fuzzyfication bloc, block ofrule 
basesanddefuzzyfication block [8], [22]. Figure 5, shows the partition of the speech universe and the 
membership functions retained for the fuzzy speed estimator (similar choices have been retained for the 
fuzzy stator resistance estimator). 


-1.00 -0.70 -0.35 0 0.35 0.70 1.00 eœ -1.00 -0.70 -0.35 0 0.35 0.70 1.00 Ae 


Figure 5. Fuzzification and partition of the speech universe for the fuzzy corrector speed estimator 
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Figure 5. Fuzzification and partition of the speech universe for the fuzzy corrector speed estimator(continue) 


For the membership functions, the triangular and trapezoidal shapes have been retained for the sake 
of simplicity. After several attempts, a partition of the speech universe with five fuzzy subsets (NB: negative 
big, PB: positive big, NS: negative small, PS: positive small, ZE: zero equal) was favored for the two 
controllers.For the treatment of the inferences, the method “Max-Min” of Mamdani was chosen for the two 
controllers, the method of defuzzification retained is that of the center of gravity.The selected rule base has 
twenty-five fuzzy rules and is described by the inference matrix shown in Table 1. According to the 
controller considered e,4e and Au represent respectively ey, (k), Ae,(k), 4@,-(k) or 


ep, (k), Aeg, (k), AR, (Kk). 


Table 1. Inference matrix defined for both fuzzy controllers 
Ae 


NB NS ZE PS PB 


NB NB NB NS PB PS 
NS ZE NS ZE PS ZE 
ZE PB PB ZE PS NB 
PS ZE PS PB NS NB 
PB PB PS NS NS NB 


5. RESULTS AND ANALYSIS 

In order to demonstrate and validate the high performance of the proposed sensorless slidingmode- 
based control system associated with FM-MRAS observer, some selected simulations are presented underthe 
software Matlab/Simulink environment. This approach is investigated by introducing a variation of load 
torque and the reference speed using a simultaneous stator resistance and rotor speed identificationof position 
sensorless induction motor to accurately estimate very low speeds. The stator resistanceis estimated 
concurrently with the speed estimation to compensate for the fluctuation in stator resistance. The study is 
applied for both induction motor operating modes; the motor mode and the regenerating mode. In this mode 
of operation, the IM issubjected to positive load torque when running at reverse speed. The parameters and 
constants of the machine are shown in Table 2. 

The accuracy of speed estimation has been verified by giving a set of steps of speed refences, two 
benchmarks are proposed and discussed in this paper as two case studies. 
- The first, designated “case 1”, shows the operation of the machine in both directions at high speeds. 
- The second, called “case 2”, discusses different operating low speed levels even in very low speed. 


Table 2. Induction motor parameters 


Designation Quantity 
Rated power 1.5kW 
Rated voltage 220V 
Rated speed 1428rpm 
Nominal frequency 50Hz 
Stator resistance 4.850 
Rotor resistance 3.8050 
Stator inductance 0.274H 
Rotor inductance 0.274H 
Mutual inductance 0.258H 
Number of pole pairs 2 

Moment of inertia 0.03 1kg/m2 
Friction coefficient 0.00114Nm.s/rd 
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5.1. Case 1: high speeds operation 

In this case, the speed curve begins with a zero speed until it varies from 0 to 20 rad/s for an interval of 
[0.75-8] s, it remains constant for an interval of [0.8-3] s. The speed then increases from 20 to 100 rad/s during 1 s, 
it remains again constant at 100 rad/s during an interval of [4]-[6] s, after that, the speed decreases to reach -100 
rad/s during 1s and remains constant until 9 s, then increases to 20 rad/s for an interval of [9-9.8] s and remains 
constant until 12 s, see Figure 6. The load torque as disturbance is varied two times from 0 to 7 Nm when the speed 
is constant for both directeanreverse sense for the two intervals respectively [4.5-5.5] s, [7.5-8.5] s, see Figure 7. 

The excellent precision of simultaneous estimation of rotor speed and stator resistance allowed to have very 
satisfactory performances in tracking even under torque steps perturbation. The zero-crossing area presents significant 
tracking and estimation errors. The stability is tested with respect to the reduction of the chattering phenomenon in 
regenerating mode of operation. In this region of operation, the motor is running in the negative speed region with 
positive load torque applied. Since the motor speed estimation depends on the online stator resistance estimation, so if 
any worsening in stator resistance estimation is fed back to the speed estimation causing instability. 

It is to notice that the estimation of the stator resistance is valid under steady state operation, see Figure 
8. As the rate of change in the stator resistance due to temperature is very low compared to the electrical time 
constants of induction motor, the FM-MRAS estimator can be disabled during transient conditions. We also 
notice that the resistance changes slightly when the load rises and the speed crosses the zero zone, and this is 
due to the increase in temperature at these moments. 


tracking speed error 


error [rad/s] 


E 2.6 
time [s] 
estimation speed error 


speed [rad/s] 


error [rad/s] 


“oO 2 4 


6 6 8 12 
time [s] time [s] 


Figure 6. IM speed tracking and estimation in different operating modes 


torque [N.m] 
curent [A] 


6 
time [s] time [s] 


Figure 7. Electromagnetic and load torques and their influence on stator currents 


stator resistance [Ohm] 
flux norm [wb] 


6 6 
time [s] time [s] 
Figure 8. Stator resistance estimate with respect to the loadtorque and norm flux 
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5.2. Case 2: very low speeds operation 

In this case we applied a benchmark of very low speed levels with testing the induction machine 
reponse in both motor and regenerating modes. The IM speed is taken as zero for the interval of [0-1.5] s, it 
then increases to 10 rad/s during 0.1 s and then remains constant until it reaches 3.5 s, the speed varies from 
10 rad/s to -10 rad/s for the interval of [3.5-3.6] s and remains constant until it reaches 4.9 s, then it increases 
from -10 rad/s to -2 rad/s during 0.1s, after that it increases again until it reaches just 1 rad/s for the interval 
of [5.3-6.1] s and keeps this value during 1.7 s then it decreases again to reache -5 rad/s at time 7.9 s and 
remains constant for the interval [7.9-9.3] s then it increases to reach the zero speed and remains constant for 
the interval of [9.4-9.9] s then it increases to 5 rad/s during 0.1 s to be constant until it attains 12 s, see Figure 
9. At the same time, the load torque is varied three times for various intervals: [2.2-3] s, [4.2-5.5] s and [8.3- 
9.2] s corresponding to variation from 0 to 3 Nm, see Figure 10. 

It should be noted that through the given tracks, we notice that the induction machine operates as a 
generator during both periods [4.2-5.5]s and [8.3-9.2]s respectively such that the IM torque and the speed 
have not from the same signs. The simulation results show an insignificant tracking and estimation errors 
even in the zone of unobservability (instability) where the speed crosses the zero zone.In order to investigate 
the performance of the controlled system for parameter variations in stator resistance. The Figure 11, presents 
various changes in stator resistance. 


tracking speed error 
T y T 


error [rad/s] 


0 2 4 8 10 1 


R 


6 

time [s] 

estimation speed error 
T f T 


speed [rad/s] 


error [rad/s] 


6 6 
time [s] time [s] 


Figure 9. IM very low speed tracking and estimation in different operating modes 
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4 
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Figure 11. Stator resistance adaptation with respect to the load torque and flux 
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At the induction machine start-up, we notice stator resistance estimation overshoots whose decrease 
when the stator resistance increases, see Figure 12. To test the robustness of tocontolled system we applied 
80% increase than 40% decrease in stator resistance respectively, see tracking and estimation speed errors 
Figures 12 and 13 (a)-(d). The simultaneous parallel estimation of the stator resistance and rotor speed 
enables the control system to eject perturbations and parameter variations. The performance results obtained 
for the proposed control strategy are acceptable and useful compared to other existing works [9], [23]-[26]. 


bd 
© 
i 


o 
o 

3 
u 


speed eror [rad/s] 


speal ator [rad/s] 
o 
spead aor [rads] 


ò 


Figure 13. These figures are;(a) tracking speed errors for +80% change in stator resistance, (b) estimation speed errors 
for +80% change in stator resistance, (c) tracking speed errors for -40% change in stator resistanceand, (d) estimation 
speed errors for -40% change in stator resistance 


6. CONCLUSION 

In this paper, fuzzy mutual MRAS simultaneous speed a stator resistance estimation has been drived for 
speed sensorless nonlinear sliding mode control of induction machine. The introducing of fuzzy mutual MRAS 
estimation in this approach achieved a high perfomance and low parameter variations sensivity given the difficulty 
comprises estimation at very low speeds where the fundamental excitation is low and the observer performance tends 
to be poor.The proposed sensorless nonlinear sliding mode control scheme showed a good behavior in the transient 
and steady states, with an excellent disturbance rejection of load torque and robustness against parameter variations. 
The simulation results demonstrate the effectiveness of proposed control over different operating conditions, a very 
good speed tracking with precise estimation in very low and zero speeds. Some perspectives of this present work can 
be about minimisation or elimination of oscillations in regenerating mode when the speed crosses the zero zone and 
also about using Dspace 1104 card to implement this approach in real-time. 
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